** Look for a unit root in real gas prices

calendar 1939 1 1

ALLOCATE 2016:1

environment echo 
source(noecho) DFUNIT.SRC
source(noecho) DOLS.SRC
source(noecho) PPUNIT.SRC

open data /mcr/home/m1res00/research/gasprice/gaspriceann.txt
   data(org=col) 1959:1 2009:1 year ngp pxg
close data
smpl 1959:1 2009:1

set rgp = ngp/pxg
set lrgp = log(rgp)
diff lrgp / dlrgp

smpl 1959:1 2009:1
@dfunit(lags=2) lrgp
@dfunit(lags=4) lrgp
@dfunit(lags=8) lrgp
@ppunit(lags=2) lrgp
@ppunit(lags=4) lrgp
@ppunit(lags=8) lrgp

@dfunit(lags=2, trend) lrgp
@dfunit(lags=4, trend) lrgp
@dfunit(lags=8, trend) lrgp
@ppunit(lags=2, trend) lrgp
@ppunit(lags=4, trend) lrgp
@ppunit(lags=8, trend) lrgp


HALT "Goodbye"
